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CHAPTER  I 


INTRODUCTION 


1. 1  Polar  Cap  Absorption  Event 

The  purpose  of  this  work  is  to  determine  and  analyze  the  electron 
density  variations  in  the  ionosphere  during  the  November  2,  1969 
Polar  Cap  Absorption  Event. 

Polar  cap  absorption  (PC A)  refers  to  an  abnormal  absorption  of 
high  frequency  radiowaves  propagating  through  a  disturbed  D  region 
of  the  ionosphere  at  polar  geomagne  •'titudes.  This  particular 
ionospheric  disturbance  has  been  attributed  to  solar  flare  emissions 
of  predominately  subrelativistic  protons  in  the  energy  range  of  5-30 
Mev  which  journey  from  the  sun  through  i.iierplanetary  space  and 
approach  the  earth.  In  their  approach  toward  the  earth,  the  earth's 
magnetic  field  deflects  and  focurcs  these  solar  particles  toward  the 
polar  regions  where  they  pei-etrate  the  earth's  upper  atmosphere 
and  ionize  the  atmospheric  constituents  in  the  altitude  range  cf  aoout 
45-75  km.  The  resultant  increase  in  the  ambient  electron  density  in 
the  D  region  results  in  the  detrimental  absorption  of  radiowaves 


which  greatly  limits  the  use  of  high  frequency  radio  and  radar 


communications  for  several  days  after  the  associated  solar  disturbance 
(flare)  has  passed. 

Research  on  PCA  phenomena  has  been  stimulated  not  only  by  the 
need  to  understand  the  deleterious  effects  on  communications,  but 
also  by  opportunities  presented  by  the  PCA  to  study  the  unfamiliar 
and  complex  chemistry  of  the  D  region  under  highly  disturbed  condi¬ 
tions.  One  of  the  opportunities  available  is  the  chance  to  obtain 
reliable  D  region  electron  density  profiles. 

During  the  PCA,  because  of  the  increase  in  the  ionization  of  the 
D  region,  the  absorption  of  high  frequency  radiowaves  propagating 
through  the  lower  ionosphere  is  quite  sensitive  to  the  D  region  electron 
concentration.  This  sensitivity  makes  possible  the  use  of  radio 
absorption  measurements  as  a  reliable  analytical  tool  for  determining 
the  D  region  electron  density.  This  ionospheric  parameter  along  with 
other  available  information  such  as  derived  electron-ion  production 
rates  during  a  PCA  event,  aids  in  the  formulation  and  testing  of 
chemical  models,  which  in  turn  lead  toward  a  better  understanding 
of  ionospheric  mechanisms;  these  include  ion-electron  formation  and 
neutralization,  ion-electron  recombination,  ion-ion  neutralization, 
and  electron  attachment  and  detachment  in  the  lower  ionosphere.  An 
excellent  review  of  the  PCA  phenomenon  was  given  by  Bailey  (1964) 
who  discussed  in  detail  most  of  the  important  related  ionospheric 
observations  and  our  understanding  of  them. 


1.  2  Operation  PC  A- 69 

Toward  resolving  some  of  the  questions  concerning  the  D  region 
processes,  the  Defense  Nuclear  Agency  (DNA)  initiated  an 

aeronomical  research  program  designated  Operation  PCA-69.  The 
program  coordinated  efforts  for  the  measurement  of  various  ionospheric 
changes  during  a  PCA  event  by  rocket,  satellite,  aircraft,  and  ground 
based  instruments.  The  objectives  of  these  measurements  were  to 
acquire  specific  data  on  the  ionospheric  changes  resulting  from  the 
sudden  influx  of  energetic  protons  into  the  atmosphere  and  the  subse¬ 
quent  effects  of  these  changes  on  high  frequency  communications. 

The  rocket  segment  of  PCA-69  consisted  of  a  series  of  scien¬ 
tifically  instrumented  rockets  launched  from  the  Churchill  Research 
Range  at  Ft.  Churchill  in  Manitoba,  Canada  during  the  November  2, 

1969  PCA  event.  Prior  to  the  solar  event,  eleven  rockets  were  fired 
to  provide  background  information  on  the  ionospheric  conditions  and 
to  provide  certification  of  the  equipment.  Thirty-six  rockets  were 
launched  to  provide  measurements  during  the  progress  of  the  PCA. 
Rocket  borne  instruments  measured  the  flux  of  energetic  protons, 
electrons,  and  alpha  particles  into  the  atmosphere;  the  radiation 
intensities  of  X-rays  and  Lyman  alpha;  and  atmospheric  effects  of 
ionization,  light  emissions,  and  heating  (Ulwick;  1971). 


1.  3  PCA-69  Areas  Rocket  Experiments 


As  part  of  the  rocket  segment  of  PCA-69,  the  Ionospheric  Research 
Laboratory  of  the  Pennsylvania  State  University  launched  nine  Areas 
rockets  designated  H-l  to  H-9  on  November  3  and  4,  1969  to  obtain 
data  on  electron  and  ion  densities  for  the  study  of  the  physical  chemistry 
of  the  D  region  during  the  PCA  sunrise- sunset  period.  Preliminary 
assessment  of  the  data  from  these  flights  w  as  given  by  Hale  and 
Seliga  (1970). 

All  nine  rockets  contained  a  parachute  borne  blunt  probe  experi¬ 
ment  to  measure  charged  particle  conductivities  which  can  be  used 
with  available  data  on  charged  particle  mobilities  and  meteorological 
conditions  to  determine  charged  particle  densities.  Hale  (1967)  and 
Hale  et  al  (1968)  described  in  detail  the  blunt  probe  theory  and 
measurements.  From  the  Areas  blunt  probe  measurements.  Hale 
et  al  (1971)  deduced  the  charge  particle  densities  during  the  November  2, 
1969  PCA  event.  These  data  are  used  later  in  this  work  to  aid  in  the 
interpretations  of  D  region  PCA  events. 

Except  for  H-9,  the  Areas  rocket  also  contained  a  CW  propagation 
experiment  which  was  to  be  operational  during  ascent  and  jettisoned  at 
apogee,  after  which  the  parachute  borne  blunt  probe  experiment 
functioned  on  descent  by  parachute.  The  propagation  experiment  was 
designed  to  determine  the  electron  density  in  the  low  D  region  of  the 
ionosphere  by  the  measurement  of  the  absorption  of  a  3.  65  MHz  CW 
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signal  transmitted  from  the  ground  to  the  rocket.  The  Ionosphere 
Research  Laboratory  conducted  the  experiments  f.nd  processed  the 
telemetry  data  to  obtain  the  signal  strength  received  by  the  rockets 
during  their  flight' (Melrose  and  Hale;  1970).  Of  the  eight  flights 
which  contained  propagation  experiments,  four  yielded  reliable  data 
from  which  electron  densities  could  be  deduced.  Pertinent  data  con¬ 
cerning  these  four  flights  are  given  in  Table  1-1. 

1.  4  Objectives 

In  this  work  the  absorption  data  obtained  from  the  propagation 
experiments  from  me  rocket  flights  of  H-2,  H-3,  H-4,  and  H-6  are 
used  to  determine  electron  density  profiles  using  numerical  full  wave 
techniques  based  on  the  generalized  magnetoionic  equations.  Compar¬ 
ison  of  the  deduced  electron  density  profiles  is  made  with  those  of 
Dean  (1971).  The  resultant  profiles  are  interpreted  in  light  of 
possible  processes  occurring  in  the  low  D  region  during  the  PCA 
event,  using  the  complementary  positive  ion  measurements  of  Hale 
et  al  (1971)  and  preliminary  calculations  of  production  rates  of 


Swider  and  Gardner  (1972). 


PERTINENT  DATA  ON  THE  FOUR  ARC'AS  FLIGHTS  FROM  WHICH  ELECTRON  DENSITIES  ARE 
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CHAPTER  II 


ARCAS  CW  PROPAGATION  DATA 

2. 1  CW  Propagation  Experiments 

The  CW  propagation  experiment  is  specifically  designed  for  use 
with  small  rockets  to  determine  electron  densities  in  the  D  region.  The 
basic  measurement  in  such  an  experiment  is  the  magnetic  and/or 
electric  field  amplitude  of  a  radiowave  transmitted  from  the  ground 
to  a  narrow  band  receiver  on  the  rocket  traversing  the  ionosphere.  A 
rocket  borne  telemetry  system  transmits  the  information  of  the 
received  signal  amplitude  to  a  ground  station  where  data  is  processed 
and  analyzed  to  determine  electron  densities  when  the  collision  fre¬ 
quency  profile  is  either  known  or  assumed  (Bowhill  and  Mechtly;  19G1). 

Hall  (1964)  and  Seliga  (1963)  performed  CW  propagation  experi¬ 
ments  to  find  D  region  electron  densities  by  measurement  of  a  specific 
amplitude  parameter  and  by  use  of  a  full  wave  solution  to  calculate  the 
parameter  based  on  several  electron  density  models.  The  electron 
density  profile  which  gave  a  field  parameter  comparable  to  the 
measured  value  was  taken  as  the  profile  representing  the  state  of  the 
ionosphere  at  the  time  of  the  flight.  The  work  by  Hall  and  Seliga 


demonstrated  the  usefulness  of  the  relatively  simple  CW  propagation 
technique  in  determining  electron  densities  in  the  D  region. 


2.  2  Details  of  the  Areas  CW  Propagation  Experiment  and  the 
Reduction  of  Data  for  Analysis 


The  Areas  CW  propagation  experiment  consisted  of  a  3.  65  MHz 
radiowave  radiated  from  the  ground  by  a  horizontal  half-wave  dipole 
to  a  rocket  as  it  passed  through  the  ionosphere.  The  transmitters 
were  located  at  the  launch  site.  In  the  rocket  a  single  ferrite  rod 
antenna, oriented  perpendicular  to  the  rocket  axis,  detected  the  magnetic 
field  strength  of  the  wave  in  a  plane  nearly  parallel  to  the  earth's 
.surface  and  served  as  a  signal  source  to  the  receiver.  Since  the 
rocket  spun  at  a  rate  of  about  10  to  20  Hz,  the  receiver  output  was 
modulated  not  only  by  the  ionospheric  effects  of  absorption  and  Faraday 
rotation,  but  also  by  the  spinning  action  of  the  rocket. 

An  example  of  the  signal  received  by  a  rocket  in  flight  is  illustrated 
in  Figure  2-1.  which  is  a  segment  of  the  processed  telemetry  signals 
from  the  H-6  flight.  The  frequency  at  which  the  nulls  appear  is  about 
twice  the  rocket  spin  frequency;  the  actual  difference  between  the  spin 
rate  and  occurrence  of  the  nulls  world  be  due  to  Faraday  rotation,  a 
rotation  of  the  plane  of  polarization  of  the  radiowave  as  it  propagates 
upward  through  the  ionosphere.  It  should  be  noted  that,  although 
planned,  Faraday  rotation  measurements  were  not  possible  because 
of  the  failure  of  the  ground  based  telemetry  receiver  to  produce  adequate 
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automatic  gain  control  signals  whose  variation  could  be  used  to 
indicate  rocket  spin  rate  as  a  reference.  The  maximum  signal  occurs 
at  the  points  where  the  orientation  of  the  ferrite  rad  antenna  couples 
most  effectively  to  the  magnetic  field  of  the  wave.  Therefore,  for  a 
rocket  having  i'.s  axis  aligned  nearly  along  the  straight  line  path  from 
the  transmitter,  the  maximum  signal  received  at  the  rocket  would 
coincide  with  the  maximum  magnetic  field  intensity  in  the  plane 
normal  to  the  direction  of  propagation;  this  is  the  case  for  the  Areas 
experiments. 

This  maximum  magnetic  iield  strength  was  selected  as  the 
measured  parameter  to  be  used  with  the  full  wave  solution  to  deduce 
the  electron  density  profiles.  This  choice  was  made  because  the 
maximum  was  a  consistent  point  that  could  easily  be  determined  from 
the  waveform  of  the  received  signal  and  because  the  maximum  was 
probably  the  most  accurate  value  detected  by  the  receiver.  In 
practice  it  was  found  that  the  relative  absorption  rates  calculated  from 
the  full  wave  solution  for  conditions  of  each  of  the  rocket  trajectories, 
using  the  same  electron  density  profile  and  accounting  for  expected 
variances  due  to  the  differences  in  the  angle  of  incidence,  were  not 
significantly  different  from  one  another.  In  addition,  the  calculated 
rates  were  practically  the  same  as  rates  obtained  from  the  maximum 
horizontal  magnetic  field  intensity  for  an  incident  wave  magnetically 
polarized  normal  to  the  plane  of  incidence.  This  latter  observation 
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considerably  simplified  the  otherwise  important  geometrical  factors 
involved  in  calculating  the  desired  field  parameter.  It  should  also  be 
noted  that  the  non-sensitivity  of  the  calculations  to  propagation  direction 
is  expected  at  Ft.  Churchill  because  of  the  propagation  conditions 
associated  with  the  near  vertical  geomagnetic  field. 

Plot  I  of  Figures  2- 2a  to  2- 2d  shows  the  maximum  relative 
magnetic  field  strength  as  received  by  each  of  the  four  rockets  passing 
through  the  lower  ionosphere.  It  should  be  noted  that  this  signal  is 
composed  of  three  components:  the  first  is  the  component  due  to  the 
ionospheric  absorption  of  the  radiowave;  the  second  is  that  due  to  the 
spreading  loss  (1/R);  and  the  third,  due  to  the  transmitting  antenna 
pattern.  For  the  analysis,  the  ionospheric  absorption  component  is 
required  and  must  be  extracted  from  Plot  I  for  each  flight. 

Using  the  radar  tracking  data  for  flights  H-2,  H-3,  H-4,  and  H-6, 
the  relative  spreading  loss  was  calculated  for  each  flight  and  added  to 
Plot  I  of  the  corresponding  figure.  Plot  II  of  Figures  2- 2a  to  2- 2d 
shows  the  maximum  magnetic  field  strength  after  spreading  loss  was 
accounted  for.  Any  compensation  of  loss  due  to  the  transmitting 
antenna  pattern  was  neglected  since  each  rocket  flew  at  approximately 
the  same  azimuth  and  elevation  angles  throughout  its  flight  in  the 
altitude  region  of  interest;  the  antenna  pattern  would  then  yield  an 
almost  constant  signal  loss  with  altitude.  Since  only  relative  absorp¬ 
tion  was  required  for  electron  density  determination,  ignoring  this 
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Figures  2-2ato  2-2d.  Relative  maximum  magnetic  field  strength 
measured  from  the  H-2,  H-3,  H-4,  and  H-6  telemetry  data.  Plot  I 
indicates  the  maximum  field  strength  as  received  by  each  rocket 
during  its  flight;  Plot  II,  the  field  strength  after  spreading  loss  is 
added  to  that  of  Plot  I. 
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constant  loss  would  not  influence  the  final  results.  Therefore,  this 
and  the  above  noted  observations  on  tne  results  of  calculated  magnetic 
field  intensities  for  different  propagation  conditions  allow  us  to  consider 
Plot  II  of  Figures  2~2a  to  2- 2d  as  the  relative  ionospheric  absorption  of 
the  maximum  horizontal  component  of  the  magnetic  field  intensity  of  a 
linearly  polarized  CW  signal  incident  on  the  ionosphere. 

Table  2-1  gives  the  absorption  rates  for  each  rocket  flight  in 
decibels  per  kilometer  derived  by  curve  fitting  the  points  given  in 
Plot  II.  These  measured  rates  which  are  independent  of  the  absolute 
field  strength  of  +he  wave  represent  the  reduced  data,  analyzed  by  the 
full  wave  technique. 
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CHAPTER  m 
ANALYSIS  OF  DATA 

3. 1  Analytical  Full  Wave  Technique 

The  analysis  of  the  data  obtained  from  the  propagation  experiments 
is  based  on  the  numerical  solution  of  the  governing  differential  equations 
of  a  plane  wave  obliquely  incident  from  free  space  on  a  horizontally 
stratified  ionosphere.  Budden  (1955)  presented  these  basic  magneto¬ 
ionic  equations  and  discussed  the  procedure  by  which  the  reflection 
coefficients  can  be  obtained  from  their  numerical  integration. 

Pitteway  (1965)  and  Seliga  (1966)  extended  the  numerical  solution  for 
use  Lx  analyzing  rocket  data.  Seiiga  (1966)  included  the  generalized 
formulation  of  the  governing  equations  which  accounts  for  the  dependence 
of  the  electron  collision  frequency  on  electron  energy  as  discovered 
by  Phelps  and  Pack  (1959)  and  analytically  treated  by  Sen  and  Wyller 
(1960).  A  computer  program  of  the  full  wave  technique  as  outlined  by 
Seliga  (1966)  was  written  for  analysis  of  the  data  obtained  from  the 
Areas  rocket  flights. 

The  program  may  be  altered  to  calculate  any  field  parameter 
desired,  but  for  the  analysis  of  the  Areas  data  the  maximum  magnetic 
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field  intensity  in  the  horizontal  plane  for  a  wave  propagating  in  the 
region  of  the  ionosphere  where  measurements  were  made  was  used. 

This  region  was  assumed  to  be  bounded  above  and  below  by  a  homog¬ 
eneous  medium.  Below  the  region  of  interest,  the  medium  is  charac¬ 
terized  bj'  a  refractive  index  having  a  real  part  within  3%  of  unity  and 
an  imaginary  part  not  greater  than  0.  0003.  For  the  purpose  of  this 
analysis,  such  a  region  was  considered  as  free  spacer  ths“?c  ly  ila  are 
consistent  with  those  considered  by  Budden  (1955). 

Above  the  region  of  interest,  the  electron  density  and  c"  !_ision 
frequency  were  assumed  to  be  constant  with  altitude.  This  assumption 
implies  that  any  reflections,  occurring  above  the  region  of  interest 
and  which  may  contribute  to  the  total  wave  field  in  the  region  of 
interest,  are  negligible.  For  the  Areas  experiments  this  assumption 
is  valid,  since  no  standing  waves  were  observed  in  the  wave  field  data 
(Figure  2- 2a  to  2- 2d),  indicating  that  the  amplitude  of  any  reflected 
wave  was  negligible  compared  to  the  upgoing  wave.  Therefore,  the 
upper  boundary  condition  of  constant  electron  density  and  collision 
frequency  does  not  affect  the  wave  solutions  in  the  region  of  interest, 
but  does  provide  an  important  artificial  physical  condition  for  numerical 
calculation  of  the  actual  wave  fields. 

These  physical  conditions  imposed  on  the  region  establish  the 
mathematical  boundary  conditions  by  which  a  numerical  solution  of  the 
differential  equations  can  uniquely  be  determined  for  an  incident  CW 


signal  of  known  polarization.  For  each  of  the:  ti  ghts,  the  incident  wave 
was  linearly  polarized  with  transmission  emanating  from  a  horizontal 
electric  dipole  aligned  in  the  East-West  direction.  Other  parameters 
required  for  the  full  wave  solution  include  a  monoenergetic  electron 
collision  frequency  profile,  the  geomagnetic  field,  the  wave  fvequ  :ncy, 
the  angle  of  incidence,  and  assumed  electron  density  profiles.  The 
fixed  parameters  for  the  Areas  experiments  are  outlined  in  Table  3-1. 

Since  the  single  frequency  CW  propagation  experiment  is  incapable 
of  determining  an  independent  measure  of  the  electron  collision  frequency, 
an  assumed  model,  shown  in  Figure  3-1,  was  required.  This  collision 
frequency  was  the  same  model  used  by  Seliga  (1968)  in  analysis  of  a 
low  frequency  propagation  experiment  at  Ft.  Churchill  in  the  summer 
of  1963.  This  model  is  representative  of  several  experimental 
measurements  of  collision  frequency.  Its  value  is  given  by 
V m  =  3.  5  x  106  exp  (-  )  sec-1 

where  z  is  the  altitude  in  kilometers. 

3.  2  Procedure  for  Determination  of  Electron  Density 

In  order  to  determine  the  electron  density  profile  an  iterative 
procedure  was  used  based  upon  absorption  rates  calculated  from  the 
full  wave  solution  using  appropriate  electron  density  profiles.  The 
procedure  is  outlined  below: 

1.  A  series  of  contour  plots  of  absorption  rate  versus  altitude 
for  a  given  electron  density  were  obtained  from  full  wave 
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Figure  .1 


solution  computations  (Fig.  3-2).  These  calculated  data 
correspond  to  the  absorption  of  a  linearly  polarized  wave 
at  3.  65  MHz,  incident  normally  on  an  ionosphere  having  a 
constant  electron  density  in  the  height  region  of  interest 
and  propagating  under  longitudinal  conditions.  These  con¬ 
ditions  approximate  those  for  each  of  the  rocket  flights  at 
Ft.  Churchill. 

2.  The  measured  absorption  rates  were  superimposed  over  the 
contour  plots  of  Figure  3-2.  The  intersection  of  these  curves 
with  the  contours  gave  the  initial  estimate  of  electron  density 
versus  altitude  for  the  respective  rocket  flights. 

3.  The  actual  conditions  of  each  flight  were  then  imposed  on  the 
propagation  conditions  and  full  wave  solutions  were  obtained 
to  compare  calculated  absorption  rates  with  the  measured 
rates.  Each  profile  was  then  altered  assuming  a  linear 
dependence  of  absorption  rate  on  electron  density  as  a 
correction.  The  full  wave  solutions  were  repeated  using 

the  corrected  electron  density  profile;  rates  were  again  com¬ 
pared,  and  the  procedure  was  repeated  until  the  rates  were 
within  three  significant  figures  of  each  other.  The  final 
profiles  were  then  presumed  to  be  the  electron  density  in  the 
ionosphere  at  the  time  oi  their  respective  flights. 


Electron /cm  3 


3.  3  Electron  Density  Results 

The  electron  density  profile  for  each  flight  deduced  by  the 
iterative  routine  described  above  is  shown  in  Figure  3-3.  The  error 
is  estimated  to  be  within  -  25%  of  the  deduced  values. 

The  results  of  the  rocket  flights  H-2,  H-3,  and  H-4  show  the 
night-to-sunrise  transition  for  November  3,  1969.  These  profiles 
show  an  increase  in  electron  density  with  decreasing  solar  zenith 
angle  most  likely  as  a  result  of  increased  photodetachment  of 
electrons  from  negative  ions.  Comparing  the  results  of  H-2  which  is 
a  nighttime  flight  with  that  of  H-3,  a  pre- sunrise  flight,  it  is  observed 
that  above  approximately  63  km  the  increase  in  the  results  of  H-3 
over  that  of  H-2  is  smaller  than  the  increase  below  this  height.  This 
observation,  assuming  little  or  no  change  in  the  production  rate, 
demonstrates  that  the  mechanism  causing  the  increase  in  electron 
density  with  decreasing  zenith  angle  becomes  more  important  with 
decreasing  altitude. 

Comparison  of  the  results  of  H-4  and  H-6  on  November  3  and  4 
respectively  shows  that,  although  the  electron  production  rate  due  to 
the  influx  of  energetic  protons  usually  decreases  with  the  progress  of 
the  PCA,  the  influence  of  solar  radiation  must  be  significant  at  low 
altitudes  in  order  to  account  for  the  nearly  factor  of  five  increase  in 
the  electron  density  between  around  55  and  63  km.  The  strong 
dependence  of  these  profiles  on  the  solar  zenith  angle  is  consistant 
with  other  observations  of  PCA's. 
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From  riometer  absorption  records  during  the  July  3-5,  1957  PCA 
Reid  and  Collins  (1959)  observed  that  this  solar  control  of  the  electron 
detachment  mechanism  was  not  predominately  dependent  on  visible 
radiation  as  evident  by  their  observations  that  the  steady  state  ionospheric 
nighttime  condition  was  reached  prior  to  true  sunset — the  time  at  which 
the  sun's  rays  are  blocked  by  the  earth.  This  observation  led  Reid  and 
Leinbach  (1961)  to  suggest  that  the  dominant  negative  ion  in  the  D  region 
during  twilight  was  not  C^,  since  visible  radiation  is  so  efficient  in 
detaching  an  electron  from  this  ion  that  if  it  were  present  in  large 
numbers,  the  detachment  process  would  be  quite  dependent  on  the 
visible  band  of  the  solar  spectrum.  The  work  of  Reid  (1961)  gave 
support  to  this  suggestion.  By  assuming  that  was  the  dominant 
negative  ion,  Reid  calculated  theoretical  absorption  curves  as  a 
function  of  solar  zenith  angle  and  compared  them  with  riometer 
absorption  measurements  of  a  PCA  during  twilight.  The  lack  of  agree¬ 
ment  between  the  calculated  and  measured  r  esults  led  him  to  conclude 
that  the  O 2  hypothesis  alone  was  inadequate  to  explain  observed 
twilight  variations  in  the  riometer  absorption.  To  explain  the  observed 
variation,  he  suggested  the  existence  of  an  unknown  ion  X  which  would 
require  ultraviolet  radiation  for  detachment.  As  yet,  the  identity  of 
X  has  not  conclusively  been  determined,  but  as  research  continues, 
the  behavior  of  this  unknov/n  ion  or  ions  should  become  better  known 
and  possibly  identified. 
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From  the  results  of  the  Areas  flights  H-2,  H-3  and  H-4,  it  appears 
that  if  photodetachment  is  responsible,  the  unknown  negative  ion 
producing  the  increase  in  electron  density  by  detachment  dominates 
below  60  km  and /  Or  the  detachment  process  is  more  efficient  below 
this  height.  Both  of  these  hypotheses  are  supported  by  the  electron 
density  results. 

The  Areas  electron  density  profiles  compare  reasonably  well  with 
the  results  of  Dean  (1971)  who  performed  similar  rocket  propagation 
experiments  during  the  same  period  of  the  PCA.  These  comparisons 
are  shown  in  Figure  3-4.  The  solar  control  which  is  evident  by  the 
increase  in  electron  density  with  decreasing  solar  zenith  angle  is  also 
noted  in  his  profiles  at  low  altitudes. 
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CHAPTER  IV 


INTERPRETATION  OF  THE 
ARCAS  ELECTRON  DENSITY  PROFILES 

4. 1  Continuity  Equation  in  the  Lower  Ionosphere 

The  D  region  of  the  ionosphere  is  a  very  difficult  region  to  under¬ 
stand,  being  complicated  by  a  large  concentration  of  various  species 
and  their  many  possible  reactions.  To  overcome  the  complexities, 
macroscopic  reactions  are  considered  which  essentially  lump  the 
numerous  reactions  of  the  different  species  into  equivalent  reactions 
involving  the  total  electron  and  positive  and  negative  ion  concentrations. 
The  continuity  equations  of  the  lower  ionosphere,  neglecting  any 
transport  processes,  are  the  following  well-known  classical  differential 
equations  (Whitten  and  Poppoff;  1971): 

W-  «  +  PpN-  +  faN-  -  AN  '  «dNN+  >'4-  » 

riN 

TjT  -  AN  -  epN_  -  paN_  -  CtiN_N+  (4.  2) 

dN. 

"dt  =  Q  ‘  adNN+  " «.iN-N+  (4.3) 

where 

_3 

N  =  total  electron  density  (cm  ) 

N_  =  total  negative  ion  density  (cm  ) 


<^*/>t%r?***«.a*iw^  ^  r^> 


N+  =  total  positive  ion  density  (cm-**) 


-3  -l 

Q  =  total  electron-ion  pair  production  rate  (cm  -s  ) 


fp 


=  photodetachment  rate  coefficient  (s  } 


-U 


-1. 


p  =  associative  detachment  rate  coefficient  (s  ) 


-1. 


A  =•  a+.;achment  rate  coefficient  (s  ) 


3-1 

=  dissociative  recombination  coefficient  (cm  -s  ) 


ct,  =  ion-ion  neutralization  (recombination)  coefficient 
1  (cm  -s  1) 


Assuming  quasi- equilibrium  conditions,  the  time  derivatives 
become  negligible,  and  the  effects  of  the  ionospheric  reaction  chemistry 
can  be  consolidated  into  an  effective  recombination  coefficient,  <X 
related  to  the  electron  and  ion  concentrations  by 

r  qi1/2 


eff’ 


N 

with 

OC 

and 

(4.  4) 
(4.5) 


*  = 


N 


N 


eP+  fa+*iN+ 


(4.6) 


The  parameters  Cf<e^  and  'X  are  characteristic  parameters  of  the 
ionosphere  at  a  given  time  and  condition.  These  parameters  can  be 
calculated  from  the  Areas  electron  density  results  and  available  data 
on  the  total  electron- ion  production  rates  and  ion  densities. 


4.  2  Calculation  of  the  Effective  Recombination  Coefficient 

Swider  and  Gardner  (1972)  calculated  preliminary  electron-ion 


pair  production  rates  from  proton  flux  measurements  at  four  different 
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times  during  the  PCA.  These  rates  are  shown  in  Figure  4-1  which 
seem  to  be  rel  atively  constant  with  height  in  the  altitude  range  of 
interest,  45-70  km. 

During  the  PCA  the  ionization  due  to  the  flux  of  solar  protons  is 
so  much  greater  than  any  present  under  normal  conditions  that  it  is 
usually  assumed  to  be  the  only  ionization  source  responsible  for  the 
production  of  electrons  (Bailey;  1S64).  Using  this  assumption. 

Figure  4-1  then  represents  the  total  electron  production  rates  in  the 
ionosphere  at  the  times  specified  by  Swider  and  Gardner. 

In  order  to  calculate  <Xe^,  the  total  production  rates  at  the  time 
of  the  four  Areas  flights  are  required.  These  rates  were  estimated 
by  plotting  the  production  rates,  deduced  by  Swider  a  id  Gardner,  versus 
time  and  then  extrapolating  between  the  measured  points.  The  results 
of  this  extrapolation  are  shown  in  Figure  4-2.  The  validity  of  the 
procedure  used  in  estimating  the  production  rata;:  during  the  Areas 
flights  depends  upon  the  assumption  that  the  production  rate  decreases 
monotonically  daring  the  course  of  the  PCA.  Such  a  decrease  is 
reasonable  in  light  of  experimental  evidence  indicating  that  the  flux  of 
protons,  with  energy  greater  than  3  MeV, which  are  primarily  responsible 
for  the  ionization  at  altitudes  between  45  and  70  km,  decreased  contin¬ 
uously  during  the  period  of  tne  Areas  flights.  The  proton  flux  measure¬ 
ments  during  this  PCA  from  the  Vela  Satellite  are  shown  in  Figure  4-3 
(Ulwick;  1971/  courtesy  of  S.  Singer,  Los  Alamos  Scientific  Laboratories). 
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Using  the  extrapolated  rates,  0(.  ^  for  each  of  the  four  Areas 
flights. was  calculated  and  appears  in  Figure  4-4.  The  estimated 
accuracy  of  cCe££  is  about  a  factor  of  two,  based  upon  the  known 

.f 

electron  density  errors  and  an  assumed  —  25%  error  in  the  calculated 
production  rates. 

The  plot  ofei,^  for  the  Areas  results  shows  that  the  electron  loss 
mechanisms  decrease  with  decreasing  zenith  angle.  Figure  4-4  also 
shows  a  comparison  of  cC^  with  the  results  of  others,  taken  from  a 
compilation  by  Poternra  et  al  (1970)  in  their  study  of  the  PCA  events 
of  1967  using  very  low  frequency  propagation  data.  It  is  noted  that 
OCgff  during  the  H-6  daytime  flight  and  the  H-2  nighttime  flight  compare 
reasonably  well  with  the  daytime  and  nighttime  results  of  Poternra 
et  al.  Between  around  58  and  63  km  the  nighttime  value  °f®^e££  was 
around  10^  -  104  greater  than  the  daytime  value. 


4.2.  1  Interpretation  of  Nighttime  Results 

Swider  etal  (1971)  interpreted  nighttime  effective  electron  loss 
coefficient  results  obtained  during  a  PCA  that  occurred  on 
November  18,  1968.  They  found  that  the  electron  density  behaves 
according  to  NCJQ/L(A),  Z  75  km,  where  L(A)  is  the  electron  loss 
rate  due  to  attachment  by  the  processes 


e  +  °2  +  o2-o2  + 


e  +  Oz  +  N^ 


02  jj.4  x  ;o‘29(^)eXp(-4^)] 

F  -31  6  -ll 

110  cm  -sec  I 
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Accordingly,  calculations  of  dL&^  using  Swider  et  al*  a  L(A)  and  the 
estimated  Q  corresponding  approximately  to  the  time  of  rocket  flight 
H-2  (Figure  4-1,  Curve  3  Nov.  69,  0005  CST)  gives  the  results  shown 
in  Figure  4-5.  This  calculated  0^e££  compares  very  well  with  the 
nighttime  values  obtained  from  0(^  =  Q/N^  and  indicates,  as  noted  by 
Swider  et  al,  that  electron  detachment  is  inoperative  at  night  and  that 
the  primary  electron  loss  mechanism  is  three -body  attachment  to 
molecular  oxygen. 

4.  3  Calculation  of  the  Negative  Ion  Density  to  Electron  Density  Ratio 
Hale  et  al  (1971)  deduced  the  positive  ion  densities  shown  in 
Figure  4-6  from  the  Areas  blunt  probe  measurements,  Comparison 
of  Hale1  s  positive  ion  densities  and  electron  densities  shows  that 
N+»  N.  Therefore,  the  neutrality  condition,  N+  =  N  +  N,  indicates 
that  the  negative  and  positive  ion  densities  are  practically  equal.  The 
parameter  ^  can  then  be  calculated  accordingly: 

A  -  It  8^-+  (4.  7) 

Results  for  those  flights  v/here  data  coincide  (H-3  and  H-4)  are  shown 
in  Figure  4-7.'  Again  the  greater  effectiveness  cf  apparent  electron 
detachment  at  lower  altitudes  with  change  in  solar  zenith  angle  is 
observed  by  noting  the  increase  in  the  difference  between  the^  of 
H-4  and  that  of  H-3  below  58  km;  above  this  height  the  /i’ s  appear 
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r*sSr.MUjrf£?L 


Altitude  (km) 


L (A),  Electron  Loss  Rote  (Ions/sec) 

Figure  4-5  a.  Electron  loss  coefficient,  L  (A),  after  Swider  et  al  (J9?1-) 


Figure  4-5  b.  Comparison  of  nighttime  ~  from  H-2  with  ~ 

N 

L(A)  given  in  Fig.  4-5  a  and  Q  given  in  Fig.  4-1. 
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to  be  nearly  equal.  These  significant  differences  in  the  s  below 
56  km  again  demonstrate  the  sensitivity  of  this  part  of  the  D  region  to 
solar  radiation  during  the  sunrise  period. 


4.4  Estimates  of  the  Ionic  Neutralization  and  the  Dissociative 
Recombination  Coefficient 


The  positive  ion  equation.  Equation  4.  3,  under  quasi-equilibrium 
conditions  may  be  written  as 


Q 


NN+  *  <*g" 


(4.8) 


At  night  and  at  low  altitudes  the  electron  recombination  process  is 
often  assumed  to  be  small,  that  is  cLrx  .  Assuming  ' 


that  this  is  valid  and  that  N  ,  then 


Q 

N  »  -Br- 

+  L 01,1 


1/2 


(4.  9) 


or 


Q 


N, 


(4.  10) 


From  the  data  of  Swider  and  Gardner  (1972)  and  Hale  et  al  (1971), 


0^  can  be  estimated  from  Equation  4.  10.  The  resultant  values  for 

-8  _7 

0 L  are  plotted  in  Figure  4-  8.  C/L  .  ranges  between  around  10  to  10 
i  i 

3  - 1 

cm  -sec  and  may  be  associated  with  two-body  ion-ion  neutralization 


involving  water  cluster  ions.  Examples  of  such  reactions  thought  to  be 
important  in  the  low  D  region  include  (Swider  and  Keneshea;  1972): 
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a;  ,lon-lon  Recombination  Coefficient  (cm3-secH ) 

Figure  4-8.  Approximate  calculation  of  the  ion- ion  neutralization 
coefficient  {ftQ/N^}  for  the  H-3,  H-4,  and  H-5 
Areas  rocket  flights. 
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H.O+  +  X  *— «»H-0  +  _ 

.3  2 

H30’r  *  HzO  +  x”—*» 2HzO  +  _ 

H  0+  *  2H  O  +  X“-*-3H  O  + _ 

.  3  2  ,  2 

H  C+  *  3H  O  +  X~— ***4H  O  + _ 

2  2  2 

H30+  •  4KzO  +  x’-**5HzO  + _ 

+  - 

H30  •  OH  +  X  — +2?!.  ^  +  . 


where.  X  represents  non-cluster  negative  ions  thought  to  be  present 

in  the  D  region,  such  as  O",  O  ,0  ,0  ,  NO  ,  NO  ,  OONO  , 

2  3  4  2  3 

CO  -  CO  .  The  rate  constants  for  these  reactions  have  not  yet  been 
3  *  4 

measured  in  the  laboratory,  but  are  expected  to  range  in  the  neighborhood 

_<*  -7  3  _} 

between  about  1  x  10  *  to  10  cm  s  with  a  rate  that  decreases  with 

increasing  hydration  {Niles;  1973). 

The  curves  show  that  Ot  reaches  a  minimum  around  48  km.  The 

i 

increase  of  o£.  with  decreasing  altitude  below  this  heignt  is  consistent 
with  expectations,  since  0(  consists  of  two  basic  neutralization 
processes:  two-body  and  three -body,  i.  e.  , 

ot,  -  *i2  +  <4-”> 


where  and  3  are  the  two-body  and  three-body  neutralization 
coefficients  respectively. 

The  variation  of  0^-  with  height  depends  on  whether  two-body  or 
three-body  recombination  dominates  and  on  the  ionic  species.  Under 


normal  atmospheric  conditions,  the  pressure  dependence  of  the  three- 


body  process  should  dominate  at  altitudes  where  the  pressure  exceeds 

1  mm  of  Hg  {Hasted;  1964).  Using  the  CIRA  1965  atmospheric  model, 

the  altitude  where  this  occurs  would  be  around  48  km.  Thus,  we 

would  expectOt  to  be  pressure  dependent  below  this  altitude 
i 

«*i  **  ^3),  increasing  at  an  exponential  rate,  probably  with  the  local 
scale  height.  Above  48  km  one  would  expect  two-body  recombination 
to  be  the  dominate  ion-ion  neutralization  process  ( oC-^-OC 2).  Further, 
if  the  ionic  constituents  remain  the  same  above  this  altitude,  we  would 
expect to  be  approximately  constant,  independent  of  pressure  and 
thus  altitude.  The  results  of  Figure  4-8  show  thatOf.  increases  with 
altitude  above  48  km.  Given  the  results  of  Narcisi  (1971)  who 
measured  positive  and  negative  ion  mass  densities  during  the  PCA, 
this  may  be  interpreted  in  several  ways;  either  ion-ion  neutralization 
is  dominated  by  a  combination  of  cluster-cluster  and  non-cluster-cluster 
neutralization  processes  that  might  change  with  height  as  the  ionic 
constituents  change,  or  Equation  4.10  is  invalid  and  electron 
dissociative  recombination  is  significant.  With  available  data  this 
first  hypothesis  cannot  be  adequately  tested.  However,  a  chemistry 
■  scheme  is  presented  in  the  next  section  to  demonstrate  a  method  of 
analysis  and  show  hew  minor  constituents  can  play  an  important  role 
in  the  actual  ion-ion  neutralization  in  the  region. 

If  48  km  is  considered  as  the  height  where  t/L  -j**  ^en  ^rorn 

Figure  4-8  and  Equation  4.  1  i  5-8x10  ^  cnA  *.  This  would 
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imply  that,  if  cluster  ions  dominate  entirely,  the  rate  for  cluster- 
cluster  reactions  is  of  this  order.  We  note  that  this  rate  has  not  yet 
been  measured  for  any  expected  reactions  in  the  laboratory  and  is  a 
consideration  which  would  require  further  study.  It  is  also  noted  that 
this  lower  estimate  for  °*i2  is  about  an  order  of  magnitude  less  than 
the  estimated  rates  used  by  Swider  and  Keneshea{1972)  in  their  compu¬ 
tations,  simulating  D  region  behavior  during  PCA  events. 

4.4.  1  Minor  Constituents  in  the  Ion- Ion  Neutralization  Process 

Although  the  calculated  C/L .  from  ionospheric  parameters  measured 
during  the  PCA  between  50-60  km  is  of  the  order  of  expected  cluster- 
non-cluster  rates,  the  height  variation  requires  explanation.  Clearly, 
changes  in  ionic  constituents  are  possible.  This  possibility  will  now 
be  considered. 

To  consider  the  role  played  by  the  various  ionic  species  in  the  ion- 
ion  neutralization  process,  the  total  negative  and  positive  ions  are  each 
decomposed  into  total  water  cluster  and  non-water  cluster  components; 
that  is 

N  =  N  +  N 

c  _  n_ 

and 

N  =  N  +  N 
+  c+  n+ 

where  N  and  N  are  the  total  negative  and  positive  water  cluster 
c.  c+ 

ions  respectively;  N  and  N  ,  the  totai  negative  ai.d  positive  non- 


water  cluster  ions  respectively.  Under  these  circumstances  the 
quasi- equilibrium  positive  ion  continuity  equation  becomes 


Q  -  06  N  N  -  0 1. 


i  c  c . 

c  c+  --  + 


,  N  N  -  Qt. 
c  n+  c_  n+ 


N  N 
:+n-  C+ 


N  N  =  0 

l  n  n, 

n_n+  -  + 


(4.  12) 


where  the  sub-scripts  of  denote  the  possible  combination  of  cluster 

and  non-cluster  ions  for  which  ion-ion  neutralization  may  occur;  for 

example,  0£.  denotes  the  coefficient  associated  with  the  neutraliza- 
c-c+ 

tion  of  a  negative  cluster  with  a  positive  cluster  ion. 

Nn  N 

Lettice  k  =  and  k .  =  ^  denote  fee  fraction  of  the  total 

&  _  M  +  - 

-  N+ 

negative  and  positive  ions  which  are  non-clusters  respectively,  we  may 
rewrite  Equation  4. 12  as 


-R  <1  -  k. 

L  c-c+ 


)  (1-KjJ  +<*.  (1-k  )  k+  + 


*c_n+ 


ot.  (l-k.-)k  + 

x  n  •  - 

C-t-  _ 


k  h  N  N  =  0 

•“-“t  +  '  J  '  + 


Hence,  we  have 

0 =  «C  (l-k  )  {l-kj +0^  (1-k  )k.  + 

c_c.  *  *  c  n,  -  + 


-  + 


Oi.  (1-k  )k  kk 

+  -  n  +  - 


1-  n 
C+  - 


(4.  13) 


-,  so  written,  reveals  the  relative  importance  that  each  of  the  ionic 
species  plays  in  the  neutralization  procesr. 
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In  particular,  if  the  non-cluster  ions  are  minor  constituents, 

-7 

i.  e.  k . ,  k_  <  0. 1  and  CL^  ft*  10  ,  then  the  last  term  in  Equation 

r  n_n+ 

4.  13  must  be  negligible,  considering  the  CL  obtained  previously  in 

i 

Figure  4-8.  Based  on  the  assumption  tha*  *:  ,  k  ^  0.  1,  then 
0.81  ^  (l-k_)(l-k+)  1, 

0  <  (l-k+)k_  0.09, 

0  <■  (l-k_)k+  <  0.08. 


From  Figure  4-7  we  may  estimate  the  magnitude  of  CC. 


“*8  '  3  —  1 

and  Of.  to  be  about  10  -  10  cm  s  ;  from  the  previous 

V- 

-9  -8  3  -1 

section  CX  was  estimated  to  be  a'jout  5x10  -  10  cm  s  . 

Vs 

Using  these  estimates,  we  can  see  that  the  first  three  terms  in 
Equation  4. 13  can  be  of  the  same  order  of  magnitude,  and  as  a  result 
the  effective  ion- ion  neutralization  would  be  dependent  or  f'e  fractional 
concentration  of  cluster  and  non-cluster  ions.  \  he.ce  resulfs  require 
further  study  and  test  as  more  specific  dais  on  cluster  ions  and  t’  exr 


chemistry  become  available. 


4.  4.  2  Electron  Dissociative  Recombination  Effects 

In  Section  4.  4  we  hypothesized  that  if  the  chemistry  is  the  same 
in  the  height  region  of  interest,  then  a  possible  explanation  for  the 
increase  in  (X.  with  height  is  that  dissociative  recombination  is 
important  and  must  be  considered  in  the  positive  ion  continuity  equation. 
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To  consider  such  a  possibility,  the  values  of  the  dissociative 

§ 

recombination  coefficient  required  to  ba'.ance  the  quasi- equilibrium 

positive  ion  equation,  using  the  Areas  electron  and  charged  particle 

densities  and  production  rates,  were  calculated  for  various  assumed 

constant  ion-ion  neutralization  coefficients.  The  curves  are  shown  in 

Figures  4-9  and  4-10  for  H-3  and  H-4  respectively. 

If  g£d  is  significant  enough  to  be  considered  in  balancing  the 

positive  ion  equation  and  if  the  chemistry  is  relatively  constant  in  the 

region  of  interest,  then  reasonable  values  for  o4»d  arid  0 C  ^  would  ba 

those  combinations  which  would  yield  constants.  Given  the  errors  in 

the  data,  such  a  region  is  shown  blocked  in  Figures  4-9  and  4-10.  For 

H-3,  the  region  is  above  about  55  km  with  04.^  and  oi^oi  the  order  of 

4-Sxlo”6  cm3  s“ 1  and  3-4xl0~8  ernes'1  respectively.  For  H-4  the  region 

is  above  54  km  with  0 C  and  oi  .  of  the  order  of  3-5x10  cm  s  and 

d  l 

1-2x10  8  cm3  s  *  respectively.  It  is  noted  that  these  values  of  are 
consistent  with  these  considered  applicable  for  this  region  for  a  complex 
water  cluster  ion  chemistry  (Aiken  et  al;  1972).  T.he  known  possible 
positive  cluster  ion-electron  dissociative  recombination  reactions  applicable 
to  this  portion  of  the  D  region  are  given  below  for  comparison  with  the 
deduced  Ot  hi  of  this  section  (Leu  et  al;  1973). 


Altitude  (krr.) 


f  igure  4-10.  Possible  combinations  of  CL d  and  Qn  necessary  to 

balance  the  positive  ion  equation  using  the  H-4  Areas 
data. 


t 

i\ 

k 


H  O  +  e-**H  +  H  O 
3  2 


(1.0*0.2x10  ^  cm^  s  *@540°K) 


H30+-  HzO  +  e  -vH  +  2H?0  (2.  2±  0.  4xl0‘6  cm3  s“ 1  @  415°X) 

H30+‘  2H20  +  e-*H  +  3H20  (3.  8  ±  0.  6xi0'6  cm3  s' 1  @  300°K) 

H30+*  3HzO  +  e  -*H  r  4H20  (4.  9-0.  8xl0'6  cm'3  s" 1  @  300°K) 

H^0+'4H20  +  e~^H  +  5H20  (6.  0  ±  1.  2xl0"6  cm3  s'1  @  205°K) 

o. 


H^O  •  5H20  +  e—*H  +  6H20  (7.  5—1.  5x10”^  cm^  s“  *  @205°K) 


H30+  •  6H20  re-+H+  7H20  (£l0-5  cm3  s"1  @  205°K) 


Other  examples  of  possibly  important,  but  still  unmeasured. 


recombination  reactions  involving  hydrated  ions  include: 


Q  *  H  O  +  e  -»-0  +  H  O 
2  2  2  2 

H  0+*  OH  +  e  — *-H  O  +  H  O 
3  2  2 

NO+  *  H  O  +  e— ►NO  +  H.,0 
2  2 


\+  . 


NO  *  2H  O  +  e  — *»NO  +  2H_0 


NO  *  3H  O  +  s— WNO  +  3H?0 

b  M 


Further  consideration  of  Figures  4-9  and  4-10  indicates  significant 


differences  foi  the  requirement  of  O C  between  H-3  and  H-4.  Below 

d 


56  km  and  for  constant  0 d.,  the  results  of  H-3  show  that 
-5.  .  ,  -4  3  _ l 

10  £  <X  <^10  cm  s  is  required  and  that  0£,  increases  with 


d  "  ’ - - d 

decreasing  height;  for  H-4,  3xl0'7jb°Cd^  6xl0"6  cm3  s"1  is  required 

and  0(_d  decreases  with  decreasing  height. 

The  large  differences  between  these  two  results  can  possibly  be 


attributed  to  the  D  region  chemistry  during  the  sunrise  period.  Reid 
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(1961)  sho wed  from  measured  riometer  observations  that  for  zenith 


o 


angles  greater  than  approximately  99  the  D  region  is  at  its  nighttime 

,o 


quasi- equilibrium  state,  and  for  angles  less  than  approximately  87 


the  D  region  is  at  its  daytime  state.  Hence,  the  results  for 


v  U 

H-3  (a»=  99.  1  )  depicts  the  nighttime  D  region  chemistry;  while  H-4 


(  87.  7°),  the  daytime  chemistry.  From  this  one  can  suggest 


.  -6  3 

that,since  rates  for  0 C  are  usually  greater  than  10  cm  s  and 

d 


since  there  is  a  tendency  for  the  rates  to  increase  with  an  increase 
in.  the  mass  of  positive  water  cluster  ion,  the  nighttime  chemistry  as 
evident  by  the  results  of  H-3  is  very  much  dominated  by  heavy  water 
cluster  ions,  while  during  the  day  the  water  cluster  ions  are  either 
lighter  in  mass  or  become  less  of  an  influence  in  the  D  region  than 
at  night.  This  suggestion  is  offered  as  an  explanation  for  the  large 


differences  between  the  s  for  "he  H-3  and  the  H-4  results 


at  the  lower  heights;  it  is  further  noted  that  such  a  suggestion  is 
open  to  question ,  but  as  further  data  become  available  we  may  be 
better  able  to  support  or  discredit  it. 
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CHAPTER  V 


CONCLUSIONS 


Electron  densities  were  deduced  from  the  H-2,  H-3,  H-4  and 
H-6  Areas  propagation  experiments  conducted  during  the  November  2, 

1969  PCA  event.  The  Areas  results  compare  reasonably  well  with 
the  electron  density  profiles  obtained  by  Dean  (1971). 

Using  available  production  rates  and  charged  particle  densities, 
the  ionospheric  parameters,  the  effective  recombination  coefficient 
and  the  negative  ion  density  to  electron  density  ratio  were  calculated. 
These  calculations  suggested  an  increasing  importance  of  electron 
detachment  during  the  night- sunrise  period  with  decreasing  height. 

The  nighttime  results  from  rocket  H-2  for  the  effective  electron  Iocs 
coefficient  compared  well  with  other  measurements  of  Swider  et  al 
(1971)  who  showed  that  electron  detachment  is  inoperative  at  night  and 
that  three-body  attachment  to  molecular  oxygen  is  most  likely  the  primary 


loss  mechanism  for  electrons  below  around  75  km. 

Estimates  of  the  ionic  recombination  coefficient,  Ot;,  were  made 


based  on  the  assumption  that  electron  recombination  was  insignificant 
in  the  altitude  range  of  about  45  to  60  km.  The  values  of  Cc  £  were 
found  to  lie  within  a  range  which  is  expected  to  apply  to  water  cluster 


ion  chemistry.  The  increase  in  o[^  below  48  km  may  be  evidence  forOd. 
being  dependent  on  pressure,  thereby  introducing  the  possibility  that 
three-body  neutralization  processes  may  begin  to  dominate  two-body 
pr ocesses below  this  height.  Possible  explanations  for  the  increase 
of  04  ^  with  height  above  48  km  were  also  giver.. 

One  explanation  was  the  possible  dominance  of  water  cluster  ions 
in  the  neutralization  processes  which  might  change  with  height.  A 
chemical  scheme  which  considered  the  possible  two-body  neutralization 
processes  between  cluster  and  non-cluster  ions  revealed  the  complex 
makeup  of  the  overall  effective  ion- ion  neutralization  coefficient  which 
can  possibly  explain  the  observed  increase  in  0^^. 

Another  explanation  was  that  dissociative  recombination  is 
important  in  the  low  D  region.  A  consideration  of  this  possibility 
revealed  that,  if  electron  recombination  is  important,  its  rate 
coefficient  would  be  of  an  order  consistent  with  reactions  involving 
cluster  ions.  Further,  it  was  suggested  that  at  night  the  water  cluster 
ions  in  the  low  D  region  would  be  heavier  than  those  involved  in  the 
daytime  chemistry. 

These  discussions  of  the  electron  density  results  were  based 
primarily  on  available  data  of  electron  production  rates  and  ion 
densities.  As  more  data  concerning  the  November  2,  1969  PCA 
becomes  available,  the  electron  densities  should  contribute  to  further 
understanding  of  D  region  chemistry  during  disturbed  polar  cap 
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conditions. 

In  addition,  the  electron  density  results  should  prove  particularly 
useful  in  testing  the  validity  of  detailed  computer  models  that  include 
the  effects  of  all  known  reactions,  production  rates,  and  transport 
effects  in  the  low  D  region. 
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